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The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured electrons 
with 0.3 < pt < 9 GeV/c at midrapidity (\y\ < 0.35) from heavy flavor (charm and bottom) decays 
in Au+Au collisions at y / ijv77 = 200 GeV. The nuclear modification factor -Raa relative to p+p 
collisions shows a strong suppression in central Au+Au collisions, indicating substantial energy loss 
of heavy quarks in the medium produced at RHIC energies. A large azimuthal anisotropy, V2, with 
respect to the reaction plane is observed for 0.5 < pt < 5 GeV/c indicating substantial heavy flavor 
elliptic flow. Both Raa and V2 show a px dependence different from those of neutral pions. A 
comparison to transport models which simultaneously describe Raa(pt) and V2(pt) suggests that 
the viscosity to entropy density ratio is close to the conjectured quantum lower bound, i.e. near a 
perfect fluid. 

PACS numbers: 25.75.Dw 



Experimental results from the Relativistic Heavy Ion 
Collider (RHIC) have established that dense partonic 
matter is formed in Au+Au collisions at RHIC 0, 0, IE 0] • 
Strong suppression observed for 7r° and other light 
hadrons at high transverse momentum (px) [E IE S Hi 
indicates partonic energy loss in the produced medium. 
The azimuthal anisotropy V2(px) [E E3] provides evi- 
dence that collective motion develops in a very early stage 
of the collision (r ^ 5 fm/c). in accordance with hydrody- 
namical calculations [Til Il3 | . The comparison of v 2 with 
several such models suggests 13J, LLJ, [l5[ that the matter 
formed at RHIC is a near-perfect fluid with viscosity to 
entropy density ratio rj/s close to the conjectured quan- 
tum lower bound [l6| . Energy loss and flow are related to 
the transport properties of the medium at temperature 
T, in particular the diffusion coefficient D oc r)/(sT). 

Further insight into properties of the medium can 
be gained from the production and propagation of par- 
ticles carrying heavy quarks (charm or bottom). A 
fixed-order-plus-next-to-leading-log (FONLL) perturba- 
tive QCD (pQCD) calculation |ljj describes the cross 
sections of heavy-flavor decay electrons in p+p collisons 
at y/s = 200 GeV within theoretical uncertainties [lij ]. 
In Au+Au collisions the total yield of such electrons 
was found to scale with the number of nucleon-nucleon 
collisions as expected for point-like processes [l9l |. En- 
ergy loss via gluon radiation is expected to be reduced 
for heavy quarks due to suppression of forward ra- 
diation, thus increasing their expected thermalization 
time [2fJ, |2l|, |22j. Consequently, a decrease of high px 



suppression and of v 2 is expected from light to charm 
to bottom quarks, with the absolute values and their px 
dependence sensitive to the properties of the medium. 
In contrast to these expectations a strong suppression of 
heavy-flavor decay electrons was discovered for 2 < pt < 
5 GeV/c [23L I24L together with nonzero electron v 2 for 
Pt < 2 GeV/c [H. 

This Letter presents px spectra and the elliptic flow 



amplitude v^ F of electrons, (e + + e~)/2, from heavy- 
flavor decays at midrapidity in Au+Au collisions at 
^snn = 200 GeV. An increase in statistics by more than 
a factor ten and reduced systematic uncertainties com- 
pared to earlier data IE IE HI greatly extend the px 
range both for the determination of the centrality depen- 
dence of Raa and for the measurement of Ug . 

The data were collected by the PHENIX detector ^ 
in the 2004 RHIC run. The minimum bias trigger and the 
collision centrality were obtained from the beam-beam 
counters (BBC) and zero degree calorimeters 0. After 
selecting good runs, data samples of 8.1 and 7.0 x 10 8 
minimum bias events in the vertex range |z v tx| < 20 cm 
are used for the spectra and v 2 analyses, respectively. 

Charged particle tracks are reconstructed with the 
two PHENIX central arm spectrometers, each cover- 
ing A<j) = 7r/2 in azimuth and \rj\ < 0.35 in pseudo- 
rapidity [2^| . Tracks are confirmed by matching showers 
in the electromagnetic calorimeter (EMCal) within 2a in 
position. Electron candidates have at least three associ- 
ated hits in the ring imaging Cerenkov detectors (RICH) 
and fulfill a shower shape cut in the EMCal, where they 
deposit an energy, E, consistent with the momentum 
(E/p — 1 > —2a). Below the Cerenkov threshold for 
pions (px < 5 GeV/c) electron mis-identification is only 
due to random coincidences between hadron tracks and 
hits in the RICH. This small background (< 20% at low 
Pt in central collisions, less towards high px and periph- 
eral events) is subtracted statistically using an event mix- 
ing technique. Requiring at least five hits in the RICH 
and tightening the shower shape cut extends the electron 
measurement to 9 GeV/c in px, with negligible hadron 
background for px < 8 GeV/c and a hadron contamina- 
tion of 20% for 8 < px < 9 GeV/c. The raw spectra 
are corrected for geometrical acceptance and reconstruc- 
tion efficiency determined by a GEANT simulation. The 
centrality dependent efficiency loss < 2% (w 23%) for 
peripheral (central) events is evaluated by reconstruct- 
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ing simulated electrons embedded into real events. 

The inclusive electron spectra consist of (1) "non- 
photonic" electrons from heavy-flavor decays, (2) "pho- 
tonic" background from Dalitz decays and photon con- 
versions (mainly in the beam pipe), and (3) "non- 
photonic" background from K — > eirv {K e ^) and dielec- 
tron decays of vector mesons. Contribution (3) is small 
(<10% for p T < 0.5 GeV/c, <2% for p T > 2 GeV/c) 
compared to (2). The heavy- flavor signal and the ratio of 
non-photonic to photonic electrons, i?NP, are determined 
via two independent and complementary methods de- 
scribed in detail in [l8| , where the identical detector con- 



figuration was used. At low pt (pt < 1-6 GeV/c), where 
the heavy-flavor signal to background ratio is small (S/B 
< 1), the "converter subtraction" method is used which 
employs a photon converter of 1.67% radiation length 
(Xq) installed around the beam pipe for part of the run. 
The converter multiplies the photonic background by a 
known, nearly px independent factor i? 7 ~ 2.3. The pho- 
tonic background can then be determined by comparing 
the inclusive electron yield with and without the con- 
verter. For higher pt, where S/B is large, the "cocktail 
subtraction" method [23| is used. Here the background 
is calculated with a Monte Carlo hadron decay genera- 
tor and subtracted from the data. At low pt the domi- 
nant background source is the 7r° Dalitz decay, which is 
calculated for each centrality using measured pion spec- 
tra [fi^l^] as input. In good agreement with measured 
data 8J], the spectral shapes of other light hadrons h (77, 
p, lo, 0, 7/) are derived from the pion spectrum assum- 
ing a universal shape in mj = yp^-\- rr? h with a fixed 
constant ratio at high pt- Photon conversions in the 
beam pipe, air and helium bags (total: 0.4%A ) are also 
included, along with background from K e ^ decays and 
both external and internal conversions of direct photons 
which are important for pt > 4 GeV/c. The agreement 
within the systematic uncertainties in the overlap region 
0.3 < pt < 4 GeV/c of these two methods demonstrates 
that the absolute value of photonic backgrounds in the 
PHENIX aperture is well-understood. 

The V2 of inclusive electrons, v™ , is measured as 
v l 2 nc = (cos[2(0 - 3>r)])/<jr [HI, where $ R is the az- 
imuthal orientation of the reaction plane measured with 
the resolution <tr using the BBC |9j. Since <tr is cen- 
trality dependent, vi is determined for narrow centrality 
bins (10%) and then averaged to calculate V2 for mini- 
mum bias events. The V2 of random hadronic background 
is subtracted statistically as described in 25 1. 



The 



non-7 



of non-photonic electrons is obtained by 



subtracting the photonic electron as: v 



non-7 



(a 



i?Np)f™ c — vJ)/Rnp- Here v\ is calculated via a Monte 
Carlo generator that includes 7r°, 77, and direct photons. 
The measured i^Cpt) of ir^,!? and [tj [2!| is used as 

input, assuming — , ^2 = v 2 > an d v^'™ ' 7 = 0- 
A direct measurement of using the converter subtrac- 
tion method confirms the calculation within statistical 



uncertainties. The resulting v™" -7 has a small contri- 
bution from K e 3 background which is simulated and sub- 
tracted to obtain uip of heavy-flavor decay electrons. 

Three independent categories of systematic uncertain- 
ties are considered. (A) The inclusive electron spectra in- 
clude uncertainties in the geometrical acceptance (5%), 
the reconstruction efficiency (3%), and the embedding 
correction (<4%). (B) Uncertainties in the converter 
subtraction are mainly given by the uncertainty in R 1 
(2.7%) and in the relative acceptance of runs with and 
without the converter being installed (1%). (C) Un- 
certainties in the cocktail subtraction rise from 8% at 
Pt = 0.3 GeV/c to 13% at 9 GeV/c, dominated by sys- 
tematic errors in the pion input and, at high pr, the 
direct photon spectrum. The V2 measurement includes a 
systematic uncertainty of 5% due to the reaction plane 
uncertainty. 

Figure [T] shows the invariant pt spectra of electrons 
from heavy-flavor decay for minimum bias events and 
in five centrality classes. The curves overlayed are the 
fit to the corresponding data from p+p collisions [Hj 
with the spectral shape taken from a FONLL calcula- 
tion [13] and scaled by the nuclear overlap integral (Taa) 
for each centrality class Q . The insert in Fig. [1] shows 
the ratio of electrons from heavy-flavor decays to back- 
ground. It increases rapidly with pt, exceeding unity for 
Pt > 1.8 GeV/c, reflecting the small amount of mate- 
rial in the detector acceptance which makes the accurate 
measurement of heavy-flavor electron spectra and v^ F 
possible. 

For all centralities, the Au+Au spectra agree well 
with the p+p reference at low pt but a suppression 
with respect to p+p develops towards high pt- This 
is quantified by the nuclear modification factor Raa — 
dN Au +Au/{(TAA}dcr p+p ), where dN Au+Au is the differ- 
ential yield in Au+Au and da p + p is the differential cross 
section in p+p in a given pt bin. For pt < 1.6 GeV/c, 
da p+p , is taken bin- by-bin from [l8j |. whereas a fit to the 
same data (curves in Fig.[T]) is used at higher pt, taking 
systematic uncertainties in da p + p and Taa into account. 

Figure shows Raa for electrons from heavy-flavor 
decays for two different pt ranges as a function of the 
number of participant nucleons, N paI t. For the inte- 
gration interval pt > 0.3 GeV/c containing more than 
half of the heavy-flavor decay electrons [T3| Raa is con- 
sistent with unity for all N paiIt in accordance with the 
binary scaling of the total heavy-flavor yield [l9j. For 
Pt > 3 GeV/c, the heavy flavor electron Raa decreases 
systematically with centrality, while larger than Raa of 
7T° withpx > 4 GeV/c p. Since above 3 GeV/c electrons 
from charm decays originate mainly from D mesons with 
Pt above 4 GeV/c this comparison indicates a smaller 
suppression of heavy-flavor mesons than observed for 
light mesons in this intermediate pt range. 

Figure [3] shows the measured Raa and of heavy- 
flavor electrons in 0-10% central and minimum bias col- 
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FIG. 1: Invariant yields of electrons from heavy-flavor decays 
for different Au+Au centrality classes and for p+p collisions, 
scaled by powers of ten for clarity. The solid lines are the re- 
sult of a FONLL calculation normalized to the p+p data [3] 
and scaled with (Taa) for each Au+Au centrality class. The 
insert shows the ratio of heavy-flavor to background electrons 
for minimum bias Au+Au collisions. Error bars (boxes) de- 
pict statistical (systematic) uncertainties. 



lisions, and our corresponding n° data jy, |2{|. The data 
indicate strong coupling of heavy quarks to the medium. 
While at low px the suppression is smaller than that of 
7T° , Raa of heavy- flavor decay electrons approaches the 
7T° value for px > 4 GeV/c although a significant con- 
tribution from bottom decays is expected at high px- 
The large v^ F indicates that the charm relaxation time 
is comparable to the short time scale of flow development 
in the produced medium. It should be noted that much 
reduced uncertainties and the extended px range of the 
present data permit the comparisons of Raa and v-i of 
the heavy and light flavors. 

More quantitative statements require theoretical guid- 
ance. Figure [3] compares the Raa and V2 of heavy- flavor 
electrons with models calculating both quantities simul- 
taneously. A perturbative QCD calculation with radia- 
tive energy loss (curves I) [30( describes the measured 
Raa reasonably well using a large transport coefficient 
q = 14 GeV 2 /fm, which also provides a consistent de- 
scription of light hadron suppression. This value of q 
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FIG. 2: Raa of heavy-flavor electrons with pt above 0.3 and 
3 GeV/c and of 7r° with pr > 4 GeV/c as function of centrality 
given by iVp ar t. Error bars (boxes) depict statistical (point- 
by-point systematic) uncertainties. The right (left) box at 
Raa = 1 shows the relative uncertainty from the p+p refer- 
ence common to all points for pt > 0.3(3) GeV /c. 
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FIG. 3: (a) 7?aa of heavy-flavor electrons in 0-10% central 
collisions compared with tt° data [f| and model calculations 
(curves I |30l ]. II [3lj . and III [32|j. The box at _Raa = 1 shows 
the uncertainty in Taa- (b) v 2 of heavy- flavor electrons in 
minimum bias collisions compared with ir° data [2S| and the 
same models. Errors are shown as in Fig. [5] 



would imply a strongly coupled medium. In this model 
the azimuthal anisotropy is only due to the path length 
dependence of energy loss, and the data clearly favor 
larger v^ F than predicted from this effect alone. 

Figure [3] also shows that the large v^ F is better repro- 
duced in Langevin-based heavy quark transport calcula- 
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tions [3lj, |32j . A calculation which includes elastic scat- 
tering mediated by resonance excitation (curves II) [3l| 
is in good agreement with both the measured Raa an d 
v 2 . This is achieved with a small heavy quark relax- 
ation time r which translates into a diffusion coefficient 
Dhq x (2ttT) = 4-6 in this model 3l|. Energy loss and 
flow are also calculated in [32] in terms of Dhq (curves 
III). While this model fails to simultaneously describe the 
measured Raa and vi with one value for Dhq, the range 
for Dhq leading to reasonable agreement with Raa or 
V2 is similar to that from [3ll |. again implying that small 
r and/or Dhq x {2-kT) are required to reproduce the 
data. Note that Dhq provides an upper bound for the 
bulk matter's diffusion coefficient D. Using the obser- 
vation [33| that D w 6 x rj/(e + p) with e + p = Ts 
at = provides an estimate for the viscosity to en- 
tropy ratio r\js w (| — 2)/4n, intriguingly close to the 
conjectured quantum lower bound 1/47T [33j]. This result 
is consistent with estimates obtained in the light quark 
sector from elliptic flow [3] and fluctuation analyses [Hf . 

The conjecture of a bound on rj/s 16| was obtained 
using the anti-de Sitter-space/conformal- field-theory cor- 
respondence [3^ . 37 1, which exploits a duality between 



strongly coupled gauge theories and semiclassical gravi- 
tational physics. Recently, such methods were ap plie d to 
estimate <?[38] and Dhq in a thermalized plasma 35^ 4(| 
41|. These authors also find a small diffusion coefficient 
D HQ x (2ttT) ~ 1. 

In conclusion, we have observed large energy loss and 
flow of heavy quarks in Au+Au collisions at ^/snn — 
200 GeV. The data provide strong evidence for the cou- 
pling of heavy quarks to the produced medium. A short 
relaxation time of heavy quarks and/or a small diffusion 
coefficient are required by the data. A model comparison 
suggests a viscosity to entropy ratio of the medium close 
to the quantum lower bound, i.e. near a perfect fluid. 
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